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ABSTRACT 

Context. In the Milky Way, most globular clusters are highly conspicuous objects that were found centuries ago. However, a few 
dozen of them are faint, sparsely populated systems that were identified largely during the second half of the past century. One of the 
faintest is ESO 37-1 (E 3) and as such it remains poorly studied, with no spectroscopic observations published so far although it was 
discovered in 1976. 

Aims. We investigate the globular cluster E 3 in an attempt to better constrain its fundamental parameters. Spectroscopy of stars in the 
field of E 3 is shown here for the first time. 

Methods. Deep, precise VI CCD photometry of E3 down to y~26 mag is presented and analysed. Low-resolution, medium signal- 
to-noise ratio spectra of nine candidate members are studied to derive radial velocity and metallicity. Proper motions from the UCAC4 
catalogue are used to explore the kinematics of the bright members of E 3. 

Results. Isochrone fitting indicates that E3 is probably very old, with an age of about 13 Gyr; its distance from the Sun is nearly 
10 kpc. It is also somewhat metal rich with [Fe/H]=-0.7. Regarding its kinematics, our tentative estimate for the proper motions is 
(Pa cosS,ps) = (-7.0+0.8, 3.5+0.3) mas yr“^ (or a tangential velocity of 382+79 km s"^ and for the radial velocity 45+5 km s"^ in 
the solar rest frame. 

Conclusions. E 3 is one of the most intriguing globular clusters in the Galaxy. Having an old age and being metal rich is clearly a 
peculiar combination, only seen in a handful of objects like the far more conspicuous NGC 104 (47 Tucanae). In addition, its low 
luminosity and sparse population make it a unique template for the study of the final evolutionary phases in the life of a star cluster. 
Unfortunately, E 3 is among the most elusive and challenging known globular clusters because field contamination severely hampers 
spectroscopic studies. 
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1. Introduction 


In the Milky Way, most globular clusters are highly conspicu¬ 
ous and rich stellar systems that were discovered centuries ago 
(see e.g. Hanes 1980; Ashman & Zepf 2008). However, a few 
dozen objects are faint and sparsely populated; they were iden¬ 
tified predominantly during the second half of the past century 
(see e.g. Harris 1976, 1996). Most globular clusters were born 
with the Galaxy, but others were captured and today are found 
associated with more or less prominent stellar streams (see e.g. 
Lynden-Bell & Lynden-Bell 1995; Ereeman & Bland-Hawthorn 
2002; Eorbes & Bridges 2010). 

One of the faintest known globular clusters is ESO 37-1 
(E 3), which was discovered by A. Lauberts in Chamaeleon in 


Send offprint requests to: R. de la Fuente Marcos, e-mail: 
rauldelafuenteinarcos@ginail. com 
^ This research note is based on observations made with the ESO 
VLT at the Paranal Observatory, under the program 078.D-0186 and 
includes data gathered with the 6.5 metre Magellan Telescopes located 
at Las Campanas Observatory, Chile (program ID CHILE-2015A-029). 
** Appendix A is available in electronic form at http://www.aanda.org. 
Tables of the individual photometric measurements are available at 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/XX/XX 


1976 (Lauberts 1976); it is still a poorly studied object and no 
spectroscopic observations have been published. Harris (2010) 
gives the position of the cluster as a(h:m:s, J2000) = 9:20:57.07, 
SC:':", J2000) = -77:16:54.8 or I = 292?270, b = -19?020, with 
My = -4.12 mag, [Ee/H] = -0.83 dex and located at a helio¬ 
centric distance of 8.1 kpc or galactocentric distance of 9.1 kpc. 
Sarajedini et al. (2007) point out that E 3 appears to be ~2 Gyr 
younger than NGC 104 (47 Tucanae), a very conspicuous 12+1 
Gyr old (McDonald et al. 2011) globular cluster. 

Here, we present and analyse deep, precise VI CCD pho¬ 
tometry of E3 down to y~26 mag. In addition, we study low- 
resolution, medium signal-to-noise ratio spectra of nine candi¬ 
date cluster members. This research note is organised as fol¬ 
lows. Section 2 is devoted to photometry. The spectroscopy is 
discussed in Section 3. The proper motions of possible bright 
members of the cluster are investigated in Section 4. In Section 
5, a polar path analysis is presented. Our results are discussed 
and our conclusions summarised in Section 6. 


2. Photometry 

On the night of 13 November 2006, E 3 was observed in the VI 
filter bands with the Very Large Telescope UT2 Kueyen and the 
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Fig. 1. y 1000-second exposure in the field of E 3. North is up, east to 


the left, and the field is 6'8 on a side. 


FORSl CCD camera. Four exposures were taken (10 s and 1000 
s in both V and I) on that night, which had nominal photomet¬ 
ric conditions and an average seeing of O'.'S, at an airmass in 
the range 1.77-1.82. The camera has a scale of 072 per pixel 
and an array of 2048x2048 pixels. An example of one of the 
6.8x6.8 arcmin^ images of E3 is shown in Fig.[T] An adjacent 
field was observed under the same conditions, offset by half a 
degree in the west direction. The standard IRApQ routines were 
applied to reduce the raw images. Using the psf-fitting routines 
of DAOPHOT and ALLSTAR (Stetson 1994) in the IRAF envi¬ 
ronment, we measured instrumental magnitudes for all the stars 
in the field. These magnitudes were converted to the standard 
system using the stars in common with Veronesi et al. (1996) 
as local standards; transformations were subsequently utilised to 
tie both the cluster and the field photometry to them (tables in¬ 
cluding the individual photometric measurements are available 
at CDS). 

The resulting colour magnitude diagrams (CMD) are shown 
in Fig. |2l for both the cluster (left panel) and the field (right 
panel). One can readily see that we can cover the cluster main 
sequence (MS) from the turn-off point (TO) down to y~25 mag, 
two magnitudes fainter than Veronesi et al. (1996). The photom¬ 
etry reaches approximately the same depth as the archival HST 
images of E 3 discussed in Lan et al. (2010). For the purpose of 
this study, we combined our photometry with that in Gratton & 
Ortolani (1987) to have a complete magnitude coverage, and to 
also sample the evolved part of the CMD (see below). Photomet¬ 
ric errors for stars with V^20 mag are ~0.03 mag (0.05 in V-I); 
for V^23 mag they amount to -0.10 mag (0.15 in V - I). 

In order to estimate star membership in the cluster, we made 
use of both the photometry of the cluster and the accompanying 
offset field. First of all, we performed a completeness analysis. 


^ IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in 
Astronomy (AURA) under a cooperative agreement with the National 
Science Foundation. 
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Fig. 2. CMD for E 3 (left panel) and the comparison field (right panel). 
Table 1. Results of the completeness test 



Cluster 

Field 

y 

comp 

comp 

(mag) 

(%) 

(%) 

20.0-23.0 

100 

100 

23.5-24.5 

90 

100 

24.5-25.0 

64 

90 

25.0-25.5 

17 

55 

>25.5 

>0 

20 

/ 

comp 

comp 

(mag) 

(%) 

(%) 

20.0-22.7 

100 

100 

22.7-24.0 

72 

100 

23.0-24.5 

24 

68 

24.5-25.0 

10 

32 


both for the cluster and for the field, which is summarised in 
Table [T] This analysis was carried out in the standard way (see 
Baume et al. 2007) by adding a number of stars per different 
magnitude bin and randomly across the field, and then comput¬ 
ing the ratio of recovered to injected stars, where the recovered 
stars result from the reduction of the added-star-image under the 
same conditions as in the original image. Looking at Table [U 
one can conclude that our photometry is complete down to y~25 
mag (7-24 mag) and that the cluster is affected by some crowd¬ 
ing. 

The completeness results have been subsequently applied to 
perform a statistical cleaning of the cluster CMD making use of 
the field, as described in Carraro & Costa (2007), with the aim of 
isolating statistically the genuine cluster population. Summaris¬ 
ing, the statistical technique we employ works as follows. We 
pick an offset field star and look into the cluster CMD for the 
closest star in colour and magnitude using a search ellipse, and 
we remove this star from the cluster star list. This procedure is 
repeated for all the stars in the offset field. We refer the reader 
to the cited paper for the full details on this procedure. In the 
following analysis we make use of this resulting clean CMD. 

Figure [3] shows the CMD for the clean sample as black cir¬ 
cles, and data from Gratton & Ortolani (1987) as red circles. Two 
isochrones with Z = 0.003, and 12 and 13 Gyr of age are used 
to perform a main-sequence fitting, giving (y - My)o = 15.07 
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Fig. 3. CMD for the clean sample (black circles). Photometry from 
Gratton & Ortolani (1987) is also shown (red circles). Two isochrones 
(12-13 Gyr, Z = 0.003) from the PARSEC evolutionary tracks (Bressan 
et al. 2012) are shown. The blue straggler VI, an SX Phoenicis variable 
star (Mochejska et al. 2000) and a probable cluster member, is shown 
as a blue circle. 

mag and E(V -1) = 0.47; the values of Ay and Aj from Schlafly 
& Finkbeiner (2011) are 0.93 mag and 0.51 mag, respectively. 
Using these values, our estimate for the cluster heliocentric dis¬ 
tance is 10.3+0.9 kpc. As the studied field is restricted to 6'8 on 
a side, we disregard estimating the cluster absolute magnitude 
and structural parameters from the photometry presented here. 

3. Spectroscopy 

Intermediate-resolution (/?j^ 7000) spectra of nine stars in the 
field of E 3 were collected at the 6.5m Baade Telescope, Las 
Campanas Observatory, on 28-29 January and 22 February 
2015. The IMACS spectrograph was used at f/4, with the 
1200+26.7 grating and O'.'75-wide slit. The resulting spectral 
range 5200-7000 A was covered by four chips, with one short 
gap between each pair of chips. Exposure times varied between 
240s and 900s, depending on target magnitude. The four spec¬ 
tral sections were reduced individually with standard IRAF rou¬ 
tines, extracted and normalised. Before merging them, the off¬ 
sets between the four independent wavelength calibration solu¬ 
tions were measured on a twilight solar spectrum and corrected. 
We measured the radial velocity (RV) of the target stars cross- 
correlating (see Tonry & Davis 1979) their spectra with synthetic 
templates drawn from the library of Coelho et al. (2005). The 
results were corrected to heliocentric velocities and for zero- 
point offsets, estimated from the position of the strong telluric 
band observed at 6850-7000 A. The errors were computed as 
the quadratic sum of the most relevant sources of uncertainty, 
namely the cross-correlation, the wavelength calibration, and 
the RV zero-point correction (for further details, see Carraro 
et al. 2007). We estimated the temperature of the targets fit- 


Table 2. Spectroscopic results. 


ID 

RA(J2000) Dec(J2000) RV 

(h:m:s) (°:':") (kms“^) 

Teff 

(K) 

V 

(mag) 

Notes 

1 

09:20:48.0 

-77:15:15 

6+9 

5500 

13.84 

field (hot) 

2 

09:20:55.9 

-77:19:01 

-23+3 

6800 

14.55 

field (hot) 

3 

09:20:41.9 

-77:18:47 

-35+3 

5750 

14.55 

field (hot) 

4 

09:20:32.1 

-77:16:25 

-21+4 

4900 

14.58 

field (met.) 

5 

09:20:47.3 

-77:18:48 

6+4 

4600 

14.74 

field (met.) 

6 

09:20:31.2 

-77:16:33 

45+3 

4950 

14.84 

Member 

15 

09:19:23.5 

-77:14:39 

44+5 

4650 

15.59 

Member 

17 

09:19:41.6 

-77:08:07 

41+4 

5600 

15.64 

field 

18 

09:19:17.2 

-77:09:40 

28+6 

5700 

16.03 

field (hot) 


ting the wings of the Hq, line with the same synthetic spectra 
(see Fuhrmann et al. 1994) following the procedure applied in 
Moni Bidin et al. (2010). We varied the gravity and metallicity 
of the template models and we thus verified that, as expected, 
the results were largely insensitive to these parameters. The un¬ 
certainty associated with the adopted minimisation is likely 
an underestimate of the true errors because it takes into account 
only random errors (spectral noise), while sources of systemat- 
ics such as continuum definition and selection of data points for 
the fit must also play a role. From the comparison of estimates 
obtained with different model parameters, normalisation, and fit 
range we estimated that the true uncertainties should likely be of 
the order of 150-200 K. The results are given in Tabled while 
in Fig.|4]we compare the position of the stars in the temperature- 
magnitude diagram with the cluster sequence. Four stars, indi¬ 
cated with “field (hot)” in the last column of the table, are too 
hot for an old stellar cluster and they will not be considered fur¬ 
ther. Unfortunately, we found no spectral feature suitable for a 
reliable gravity estimate at our low resolution. The degeneracy 
between gravity and metallicity in the fit of the spectra could 
therefore not be fully broken and we derived the best-fit metal¬ 
licity when log g was varied in the range 0-4 dex. The metal¬ 
licity of star 4 resulted [Fe/H]> -0.35 for any value of logg, 
and the best-fit solution of star 5 indicated that it is most likely 
a dwarf star (logg > 4). Their cluster membership could thus 
be excluded and they are flagged with “field (met.)” in Tabled 
The metallicity of star 17 was also too high ([Fe/H]j^0.0) to be 
of interest here. Object 6 is likely a red giant (log^ 1.8) with 
[Fe/H]=-0.65 + 0.15, while the spectrum of target 15 is compati¬ 
ble with either a metal-poor giant (log g ~ 2, [Fe/H]= -0.7+0.2), 
or a metal-rich dwarf. 

In conclusion, only two stars in the observed sample (6 and 
15) have both temperature and metallicity fully compatible with 
cluster membership, and their RVs differ by only 1 km s“^. Tar¬ 
get 17 also has similar RV, but - although its higher temperature 
is still compatible with a horizontal branch cluster star - its high 
metallicity points to a field star. Therefore, our estimated value 
for the radial velocity of E3 is 45+5 km s“^ in the solar rest 
frame. 


4. Proper motions 

The two most robust candidates for membership in E3 have 
proper motions in UCAC4 (Zacharias et al. 2013). Star 
6 is UCAC4 064-008502 with (//^ cos S,/is) = (-6.4+1.9, 
3.3+1.7) mas yr-^. Star 15 (UCAC4064-008471) has consis¬ 
tent proper motions of (-7.6+2.0, 3.7+2.3) mas yr“^. Star 17 
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Fig. 4. Comparison between the position of the stars studied spectro¬ 
scopically (see Table O in the temperature-magnitude diagram with an 
assumed cluster sequence, Z = 0.003 and age of 13 Gyr, from the PAR¬ 
SEC evolutionary tracks (Bressan et al. 2012). 


(UCAC4 065-008924) has somewhat inconsistent proper mo¬ 
tions, (-1.6+3.5, 1.1 ±5.4) mas yr“^ Assuming that stars 6 and 
15 are part of the cluster, our tentative estimate for the proper 
motions is (pa cosS,ps) = (-7.0+0.8, 3.5+0.3) mas yr“^. A tan¬ 
gential velocity of 382+79 km s“^ is derived from the inferred 
cluster proper motion. 

The globular cluster E 3 is a low-luminosity stellar system 
with very few bright stars in its central regions. The closest star 
to the assumed centre of the cluster (VizieR, J2000 09:20:57.1 
h:m:s, -77:16:55 °:':") with proper motions in UCAC4 is 
40" from the centre. The closest putative member could be 
UCAC4 064-008511 with (-7.0+4.3, 4.1 + 13.9) mas yr'^ 63" 
from the centre. Stars 6 and 15 are located 1'4722 and 5'6358, re¬ 
spectively, from the centre. In contrast, star 17 is located 9'0846 
from the cluster centre. Van den Bergh et al. (1980) estimated 
the core and tidal radii of E 3 at 4.4 pc (1'51) and 24.7 pc (8'49), 
respectively. Being extra-tidal is another argument against star 
17 having membership in E 3. 

Eigure [5] shows the 2MASS near-infrared CMD of selected 
stars. The objects plotted have proper motions compatible with 
those of stars 6 (064-008502) and 15 (064-008471) and are lo¬ 
cated within the cluster tidal radius. The sources VI and V2 are 
putative cluster members from Mochejska et al. (2000), but V2 
has no proper motions in UCAC4 and V1 has rather incompati¬ 
ble proper motions, (0.7+5.7, 22.2+5.7) mas yr“^ 


5. Polar paths 

Most globular clusters were born with the Galaxy, but some oth¬ 
ers were accreted. Each Galactic globular cluster and satellite 
dwarf galaxy orbits the Galaxy in a plane including the object’s 
present position and the Galactic centre. This reference plane 
can be defined as the instantaneous or osculating orbital plane 
of the object. Eollowing Lynden-Bell & Lynden-Bell (1995) we 
consider the geometric loci of all possible normals to the galac- 
tocentric radius-vector for a particular object (see Appendix A 



Fig. 5. 2MASS near-infrared CMD of stars within 8:5 from the cluster 
centre and with proper motions within 2cr of our estimate for the cluster 
(black circles, labels are the UCAC4 designations). Variable stars VI 
and V2 from Mochejska et al. (2000) are also indicated (empty circles). 

for mathematical details). This defines the set of possible orbital 
poles for that object, globular cluster, or dwarf galaxy. 

For E3, the usual heliocentric galactic coordinates (l,b) = 
(292?270, -19?020) translate into galactocentric gal actic coordi- 
nates (lg,bg) = (232?550, -16?474). We solved Eq. lA.5l using a 
Monte Carlo approach and evaluated the statistical significance 
of all the crossing points of polar path intersections with that 
of E3 using the Kulldorff statistics (Kulldorff 1997). All studies 
coincide in classifying E 3 as a metal-rich globular cluster. For 
this globular cluster subset the most relevant intersection (180 (t) 
appears to be located at X = -0.42, Y = 0.50 or (/, b) = (310°, 
+35°), see Fig.[6j In these diagrams a separation of 0.1 is equiv¬ 
alent to ~ 6°. The poles in Fig.[6]are all within 2° of the average 
pole. Prospective members of this dynamical family of objects 
include, besides E 3, the second most massive Galactic globular 
cluster NGC 104 (47 Tucanae), NGC 6171 (M 107), and NGC 
6362. Additional members may be Terzan 2 and NGC 6652. 
All of them have ages in the range 12.8+0.6 Gyr and metallic- 
ity [Fe/H] in the range (-0.69, -1.02). If we compare the polar 
paths of this group of objects with those of objects commonly 
associated with the Sagittarius stream (see Figs.l6land lA.ll) . the 
dispersion of the poles is even smaller in this case. 

6. Discussion and conciusions 

Clearly, more observations are required to investigate this inter¬ 
esting cluster further. In particular, detailed spectroscopy of the 
brightest stars can provide important information on the actual 
metallicity of the cluster and its radial velocity. However, our 
preliminary spectroscopic study has revealed that this cluster is 
profusely contaminated with field stars. On the other hand, its 
low luminosity and sparse population make it a unique template 
for the study of the final evolutionary phases in the life of a star 
cluster. An extensive spectroscopic study can be used to measure 
the radial velocity dispersion that in turn can be utilised to con¬ 
strain the dynamical state of this object. In a cluster with a high 
binary fraction like E 3 (e.g. Veronesi et al. 1996), binary orbital 
motions can mask small velocity dispersions. Multi-epoch ob¬ 
servations can be used to separate single from binary systems. 
Unfortunately, E3 is among the most elusive and challenging 


Article number, page 4of|7] 







de la Fuente Marcos et al.: The globular cluster ESO 37-1 (E 3) (RN) 


known globular clusters because field contamination severely 
hampers spectroscopic studies. 

The globular cluster E3 does not have a clear horizontal 
branch and that makes its distance and age determination rather 
uncertain. Our values differ significantly from those in Sarajedini 
et al. (2007), namely 8.1 kpc and 10 Gyr, respectively. The value 
of Ay used in this study is nearly 17% lower than those used in 
previous works. Our isochrone fitting indicates that E 3 is prob¬ 
ably very old, with an age of about 13 Gyr and its heliocentric 
distance is nearly 10 kpc. We have confirmed that it is also some¬ 
what metal rich with [Ee/H]=-0.7. Eor its kinematics, our tenta¬ 
tive estimate for the proper motion is (fia cos S, fis) = (-7.0+0.8, 
3.5+0.3) mas yr“^ A tangential velocity of 382+79 km s“^ is 
derived from our determination of the cluster proper motion and 
a radial velocity of 45+5 km s"^ is obtained from our spec¬ 
troscopy, in the solar rest frame. The E 3 stellar system is one 
of the most intriguing globular clusters in the Galaxy. Its old 
age and high metallicity is a peculiar combination, only seen 
in a handful of objects like the far more conspicuous NGC 104 
(47 Tucanae) that may share a stellar stream with it. 
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Appendix A: Polar path technique analysis 

Each Galactic globular cluster and satellite dwarf galaxy orbits 
the Galaxy in a plane including the object’s present position and 
the Galactic centre. This reference plane can be defined as the 
instantaneous orbital plane of the object. Following Lynden-Bell 
& Lynden-Bell (1995) we consider the geometric loci of all pos¬ 
sible normals to the galactocentric radius-vector for a particular 
object. This defines the set of possible orbital poles for that ob¬ 
ject, globular cluster, or dwarf galaxy. We consider a system of 
galactic coordinates on the unit sphere = 1) cen¬ 

tred on the Galactic centre; the coordinates of any object in that 
system are given by 

(v, y, z) = (cos Ig cos bg, sin Ig cos bg, sin bg) , (A.l) 

where Ig and bg are the galactocentric (not the regular heliocen¬ 
tric) galactic longitude and latitude, respectively. The galacto¬ 
centric radius-vectors for an object and its pole are perpendicu¬ 
lar, therefore 

V cos Ig cos bg y sin Ig cos bg z sin bg = 0. (A.2) 

Finding z(x, y) from Eq. IA.2[ replacing it on the expression of 
the unit sphere, and reorganising we obtain 



If we now consider the Lambert azimuthal equal-area projection 
given by the expressions 


{X, Y) = {x/ vm, y/ VTTI), 

with z = 1 - - Y^, Ea. IA.BI can be rewritten as 


(A.4) 



= 1. (A.5) 


This expression provides the positions of all the possible poles 
associated with the object. The intersection (if any) of polar 
paths gives the putative pole of a group of (perhaps) dynam¬ 
ically related objects. Once the pole has been found, its po¬ 
sition on the sky relative to the Galactic centre is given by 
(Ip, bp) = (arctan(Tp/Xp), arcsin(l - Xp - Yp), see Eqs. lA.ll 
and lA.4l Polar paths of objects possibly associated with E3 are 
shown in Fig. [b] For comparison, those of objects widely re¬ 
garded as part of the Sagittarius stream, streams 8a and 8b in 
Tabl e 2 of Lynden-Bell & Lynden-Bell (1995), are plotted in Fig. 
lATl 
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Fig. A.l. (Left panel) Polar paths of objects believed to be part of the Sagittarius stream, streams 8a and 8b in Table 2 of Lynden-Bell & Lynden- 
Bell (1995). (Right panel) Close-up. 




Fig. 6. (Left panel) Polar paths of objects believed to be part of the E 3 stream. (Right panel) Close-up. 



















